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Abstract —Minimizing the cost of deploying and operating 
a Wireless Sensor Network (WSN) Involves deciding how to 
partition a budget between competing expenses such as node 
hardware, energy, and labor. Most commercial network operators 
account for Interest rates in their budgeting exercises, providing 
a financial incentive to defer some costs until a later time. In 
this paper, we propose a net present cost (NPC) model for WSN 
capital and operating expenses that accounts for Interest rates. 
Our model optimizes the number, size, and spacing between 
expenditures in order to minimize the NPC required for the 
network to achieve a desired operational lifetime. In general 
this optimization problem is non-convex, but if the spacing 
between expenditures is linearly proportional to the size of the 
expenditures, and the number of maintenance cycles is known in 
advance, the problem becomes convex and can be solved to global 
optimality. If non-deferrable recurring costs are low, then evenly 
spacing the expenditures can provide near-optimal results. With 
the provided models and methods, network operators can now 
derive a payment schedule to minimize NPC while accounting for 
various operational parameters. The numerical examples show 
substantial cost benefits under practical assumptions. 

Index Terms —Wireless sensor network (WSN), net present 
cost (NPC), net present value (NPV), cost, budget, lifetime, 
deployment. 


I. Introduction 

W IRELESS Sensor Networks (WSNs) are groups of 
nodes that collaboratively collect information on an 
area of interest. Their ability to reduce costs and save human 
lives by autonomously monitoring remote and potentially haz¬ 
ardous regions has made them an active area of research, with 
applications in smart agriculture, environmental monitoring, 
detecting faults in systems and structures, disaster monitoring, 
and battlefield surveillance [1]. Nodes consist of sensors and 
transceivers to gather data on their immediate surroundings 
and forward this data over an ad-hoc network structure to 
predefined locations for further processing. In order to min¬ 
imize the cost of covering an area of interest, nodes are 
designed with inexpensive hardware, implement low-power 
protocols such as ZigBee [2], and may use scheduling [3] 
and energy-minimizing routing [4]. By utilizing robots [5] and 
unmanned aerial vehicles [6] to replenish energy and replace 
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damaged nodes, a network operator may reduce labor costs 
while extending WSN lifetime. 

A network operator allocates a limited budget to numerous 
tasks related to building and maintaining a sensor network: 
hardware must be purchased and deployed, batteries may re¬ 
quire periodic replacement or recharging, and damaged nodes 
may need to be replaced. Allocating additional money to one 
part of the budget reduces available funds for the other parts; 
for example, adding nodes to a network increases the portion 
of the budget dedicated to node hardware, but decreases the 
money available for energy and labor. As discussed in Sec. II, 
a large body of research focuses on minimizing the individual 
WSN costs, such as focusing solely on node hardware costs 
or node energy costs, but only a few papers study how to 
minimize the overall cost when multiple different expenses 
are combined together. 

Before undertaking a project or submitting a bid, compa¬ 
nies often estimate that project’s initial investment, known 
as Capital Expenditures (CAPEX), and recurring expendi¬ 
tures, called Operational Expenditures (OPEX). Eor a WSN, 
CAPEX includes the costs of node hardware, the initial node 
energy supplies, and the labor required to initially deploy 
the network. OPEX includes the cost of replacement node 
hardware, replacement batteries, and the labor required to 
perform maintenance on the network. 

To improve the cost of a WSN we propose a framework for 
minimizing its Net Present Cost (NPC). NPC is similar to net 
present value [7], except all cash flows are considered outflows 
instead of being either inflows or outflows. NPC combines 
CAPEX and OPEX into a single cost by taking interest rates 
into account. Purchases made in the future cost less, in terms 
of the present currency value, as the operator earns interest 
by collecting revenue from the network it has built and by 
investing money elsewhere. Our framework could be used, for 
example, by a network operator that wants to take advantage 
of interest rates to reduce the cost of energy in the future. 
This requires spending more on labor in the future to visit the 
network and deliver this less expensive energy to the nodes. 
Our NPC minimization framework would, in this case, find 
the optimal balance between the money saved on energy and 
the cost of labor required to deliver it to the network. 

By deferring costs to take advantage of interest rates, NPC 
minimization can significantly reduce the total cost of a WSN 
compared to paying for all costs up-front. Our framework 
produces a schedule of expenditures that minimizes NPC; this 
schedule can be used by a network operator when estimating 
the budget of a WSN. General rules-of-thumb can be applied 
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in certain scenarios to produce near-optimal budgets. When 
non-deferrable recurring costs such as labor costs are low, the 
sensitivity of NPC minimization to the number of maintenance 
visits performed is also low, meaning that performing main¬ 
tenance as often as possible produces a near-optimal NPC. In 
addition, we show that when non-deferrable recurring costs 
are low, evenly spacing the maintenance visits can provide a 
near-optimal NPC. 

We dehne a visit as a time point where the network operator 
visits the network to perform scheduled maintenance, such 
as restoring energy to nodes in the network. Each visit has 
an expenditure associated with it, referred to as the visit 
expenditure, while the time until the next visit is called the 
visit lifetime. The visit made to initially deploy the network 
is a CAPEX expenditure, while visits made after deployment 
to restore node energy are OPEX expenditures. The network 
operator can adjust the number of visits, as well as each 
visit expenditure and visit lifetime in order to minimize 
the NPC. These parameters are interrelated, so adjustments 
are not always straightforward; for example, increasing visit 
expenditures may increase visit lifetimes and reduce the total 
number of visits required to achieve a desired operational 
lifetime. This action is only worthwhile if the reduced number 
of visits compensates for higher cost per visit. 

We propose a two-layer optimization framework for de¬ 
termining the number of visits, the visit expenditures, and 
the visit lifetimes required to minimize the NPC. The hrst 
layer of the framework is a non-convex optimization problem 
that maximizes visit lifetime when given a visit expenditure 
that is known in advance. The visit lifetime depends on 
the visit expenditure; the more money spent on energy, the 
longer the visit lifetime. Maximizing visit lifetime minimizes 
future costs by taking full advantage of interest rates. This 
optimization problem is used to derive a lifetime function for 
each visit. The lifetime function represents the relationship 
between visit expenditure and the maximum visit lifetime that 
can be achieved with that visit expenditure. Each visit may 
have its own unique lifetime function to account for changes 
in node hardware, energy, and labor costs between visits. 

The second layer of the framework optimizes visit expen¬ 
ditures and the number of visits to minimize the NPC. It uses 
the lifetime function from the hrst layer of the optimization 
framework to calculate the optimal visit lifetime for each visit 
expenditure. This ensures that the visit expenditures found by 
the second layer provide the optimal visit lifetimes. We show 
that the optimization problem in the second layer is non-linear 
and non-convex, making it difficult to hnd a globally optimum 
solution. 

To reduce the complexity of this non-convex problem, we 
show that if the total number of visits is hxed to K, and the 
lifetime function for each visit is piecewise linear (equal to 
zero until a certain visit expenditure and increasing linearly 
with slope m afterwards), then NPC minimization will be 
convex and the solution found will be globally optimal. The 
slope 771 will be equal for all visits, but the point where the 
function transitions from zero to an increasing linear function 
may differ. Assuming that a maximum of K visits are allowed, 
we provide a 0{K^) algorithm for hnding the number of visits 


that minimize the NPC. The lifetime functions are piecewise 
linear under a network model that assumes optimum data 
flows between nodes, optimum energy consumption, and that 
1-connectivity is adequate. 

Points in time where the network operator performs un¬ 
scheduled maintenance on the network, due to unexpected 
events such as hardware faults or environmental damage, 
are considered when calculating the NPC. As such events 
may occur at any given moment, we do not know when an 
unscheduled repair will occur, so we approximate the NPC 
of unscheduled payments by assuming that failures occur 
periodically, with the length of each period equal to the Mean 
Time Between Eailures (MTBE) [8] of the network nodes. The 
NPC of unscheduled payments can be reduced by improving 
the reliability of the network. Doing so increases the cost of 
node hardware, either by purchasing more robust hardware or 
by adding redundant nodes in stand-by to node locations, the 
locations of active nodes in the network. To balance cost and 
robustness, we provide a method of minimizing node hardware 
costs and the approximate NPC of unscheduled repairs. The 
method assumes every node has identical hardware and every 
node location is given the same number of redundant nodes in 
stand-by, and hnds the best type of node hardware and level 
of redundancy out of a set of possible choices. 

A brief survey of the literature relevant to minimizing 
deployment and operational costs of WSNs is provided in 
Sec. II, while our network and cost models are given in 
Sec. III. Sec. IV discusses how to maximize visit lifetime 
for a given visit expenditure by optimally dividing the visit 
expenditure between node hardware, energy, and labor costs. 
Sec. IV also demonstrates that the lifetime function, the 
relationship between visit expenditure and visit lifetime, is 
piecewise linear under our network and cost models. A general 
NPC minimization framework is proposed in Sec. V, and the 
linear lifetime function found in Sec. IV is exploited to show 
that equally spacing visits can be a near-optimal technique 
for minimizing the NPC under our network model. Sec. VI 
provides numerical results for the NPC minimization problem 
under practical assumptions. It discusses the effect of network 
lifetime, the interest rate, the MTBE of the network nodes, 
and the costs of node hardware, energy, and labor on the NPC 
of the network. 

II. Related Work 

When minimizing the overall WSN cost, NPC minimization 
accounts for the costs of node hardware, energy, and labor. 
The majority of papers available in the literature tend to focus 
on only a single type of cost. Both [9] and [10] are survey 
papers that contain techniques for minimizing the node count, 
and therefore node hardware costs, when connecting disjoint 
networks. The survey papers [11] and [12] discuss techniques 
for reducing network energy consumption and therefore en¬ 
ergy costs. Labor costs can be reduced by scheduling node 
replacement and recharging activities [13] and automating 
maintenance processes via robots [5] and unmanned aerial 
vehicles [6]. While the techniques presented in these papers 
can reduce network cost, they are not useful for operators 
seeking to allocate a budget across multiple costs. 
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There is a significant body of work on optimizing the 
deployment and operational costs of cellular networks [14], 
[15], [16], [17]. The fundamental difference with respect to 
this paper is that, once deployed, the infrastructure of a 
cellular network is static while the topology of a WSN varies 
in time. This variation results from failed sensor nodes due 
to depleted batteries and the dynamics of ad-hoc routing. 
Compared to cellular networks, WSNs require a substantially 
different model to compute the deployment of energy (i.e. 
batteries) and maintenance cycles. 

To the best of our knowledge, the only papers in the 
literature accounting for multiple types of costs in a WSN 
are the works by Misra et al. [18] and Dutta et al. [19]. 
Compared to these papers, the NPC minimization approach in 
this paper is unique in that it determines the visit expenditure, 
visit lifetime, and number of visits while taking advantage 
of interest rates. The other research that combines multiple 
types of network costs together either ignore savings from 
interest rates or fail to optimize the spacing between visits, 
both of which are critically important components of creating 
a network budget. 

Misra et al. [18] provide a method of balancing the cost 
of WSN maintenance with network performance loss. Per¬ 
formance refers to the probability of the WSN detecting an 
event; this probability, and therefore performance, decreases 
as the number of failed nodes increases. This work uses a 
probabilistic model to predict the long-term cost of a network, 
and provides an algorithm for determining the optimal trade¬ 
off between node replacement cost and performance loss. 
Their technique focuses solely on OPEX, ignoring the CAPEX 
required to deploy the network, and does not factor in savings 
due to interest rates, meaning that reductions in maintenance 
costs from returns on investments or revenues generated by 
the network are not considered when determining the node 
replacement policy. 

Dutta et al. [19] provide a strategy that considers interest 
rates when determining when to replace nodes in a network 
to minimize maintenance costs. This research determines the 
optimal year in which to replace nodes and assumes that 
maintenance costs are known in advance. It does not find 
the optimal size of each visit expenditure nor the optimal 
spacing between visits, limiting its applicability to minimizing 
a WSN’s budget. By determining the amount to spend on each 
visit, the number of visits, and the spacing between them to 
minimize NPC, we allow network operators to not only reduce 
the total WSN cost, but to also estimate the budget and cash 
flows of a WSN. 


III. Network and Cost Models 
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Fig. 1. Payments made over the course of a network’s lifetime. Scheduled 
payments pi ... pn are known in advance, while unscheduled payments 
(Ti . .. ap occur at times t\ .. .tp. 


A. Network Model 

Every network has a set of sensor node locations T, a 
set of sink locations S, and a set of potential or candidate 
relay locations M. The union of these sets is referred to as 
the universal set = TU 5 U AA. Relay nodes forward data 
from sensor nodes to the sink nodes, possibly over multiple 
hops. Not every candidate relay location will have a relay node 
placed on it; depending on factors such as the transmit range 
of other nodes and the network topology, certain candidate 
relay locations may be chosen over others. Sensor nodes, in 
addition to performing the same forwarding duties as relay 
nodes, collect and forward data from their on-board sensors. 
Sink nodes collect the data gathered by every sensor node in 
the network and store or process it. Sensor node and sink node 
locations are assumed to always have sensor and sink nodes 
respectively placed on them. 

Sink, sensor, and relay nodes placed on their respective 
locations consume energy and handle data while performing 
their roles. Each sensor node i in Tgenerates data at a rate of 
gi bits/s and consumes energy at a rate of Si J/bit when sensing 
data. Each sensor or relay on location i in AT U T consumes 
Ti{d)}Po\\. when transmitting to a location that is distance d 
meters away, and consumes Qi J/bit while receiving data. An 
edge exists from node i to another node j if the distance 
between the nodes dij is less than the maximum transmit 
distance in meters. 

When determining how many nodes to place, where to 
place them, and how much energy to allocate to the nodes 
in a network, the operator has specific goals in mind; here 
we assume the goals are to achieve an operational lifetime 
of exactly L years and 1-connectivity at minimum cost. 
While other factors such as latency and quality of service are 
important in a number of situations, we assume real-time data 
acquisition is not critical and that the network will be lightly 
loaded, so that such factors are not a priority. 


In this section we provide a network model and a cost 
model for the NPC minimization problem. Sec. III-A describes 
the types of locations present in the area of interest, their 
properties, and how they relate to one another to form a WSN. 
In Sec. III-B we present a model for NPC that takes into 
consideration the number of visits, the visit expenditures, and 
the visit lifetimes. 


B. Cost Model 

Pig. 1 illustrates our cost model. We consider two types of 
payments: scheduled payments for performing routine tasks 
such as restoring energy to nodes, and unscheduled payments 
for reacting to unexpected issues such as hardware faults. 

Scheduled payments are divided into events called visits. 
Assuming that a total number of K visits are made, there 















4 


is an initial visit expenditure pi that accounts for initially 
deploying the network, and a number of subsequent visit 
expenditures p 2 ■ ■ ■ Pk for maintenance. The time between pk 
and Pk+i is the visit lifetime Ik', as will be discussed in Sec. IV, 
the expenditure pk is optimally divided between hardware, 
energy, and labor costs in such a way as to maximize Ik- 
The operational lifetime of the network is L = h- When 
minimizing the NPC of a network, one has to determine each 
visit expenditure pk, the total number of visits K, and the visit 
lifetimes Ik- 

Unscheduled payments occur to handle network failures 
that cannot be prevented through routine maintenance, such 
as those caused by hardware faults, environmental hazards, or 
accidental damage. We assume that F failures occur at times 
ti - - - tp, with failures costing ai - - - ap financial units each 
to repair. 

Minimizing the cost of a network with a long operational 
lifetime entails determining whether money is best spent 
immediately, or later after it earns interest at a rate of u > 0. 
In other words, one way to minimize costs is to minimize the 
Net Present Cost (NPC) of the network. Net present cost can 
be expressed as 


NPC 


Pi + 111=1 
Pi + Y.k=2 


(l+T))* 


£fc 

ir 


(1+Tl)^n = l 


+ 

^ Z^n=l 


(l+i;)* 


K = 1 
K >2- 

( 1 ) 


We assume that a relationship exists between the visit 
expenditure pk during visit k, and the visit lifetime Ik- In 
other words 


h = fk{pk)- 


( 2 ) 


Each lifetime function fk{pk) may be unique for each visit k- 
As we do not know when failures occur, we propose a 
method to approximate failure times for the purpose of this 
analysis. We assume that the network’s Mean Time Between 
Failures (MTBF) uj is known, that uj does not vary over the 
network’s operational lifetime L, and that the cost a of failures 
does not change. As discussed at the end of Sec. III-A, we 
assume that 1-connectivity is adequate for the network. In 
this case, a single node failure potentially disconnects the 
network, meaning that w is the same as the MTBF of the 
network nodes. With our approximation, = nui and (T„ = cr, 
Vn £ [1, [L/wJ], meaning that the NPC can now be written 
as 


NPC = 


Pi 


(j 

m—1 (1+1))”'^ 

Pk 

K = 1 

>LL/ci;j (J 

•n—1 

K>2 


(3) 


It is important to point out that minimizing overall NPC 
requires two layers of optimization. In Sec. IV, a first op¬ 
timization is used to determine how an expenditure pk is 
optimally divided between hardware, energy, and labor costs 
in order to maximize Ik- The results of this optimization 
establishes the function Ik = fk{Pk)- In Sec. V we then use 
the function fk{Pk) in n second optimization that minimizes 
overall NPC by optimizing the visit expenditures pk- 


IV. Framework for Optimizing Visit Fifetime 

From (1) we see that minimizing NPC involves maximizing 
visit lifetime, Ik- We also assume in (2) that a relationship 
Ik = fk{Pk) exists. In this section we provide an optimization 
problem that maximizes Ik for a given pk by optimizing the 
information flow rates in the network and the division of pk 
between hardware, energy, and labor costs. This optimization 
is formulated as a non-linear, non-convex, yet continuous 
problem in Sec. IV-A, and is turned into a non-continuous 
Mixed Integer Program (MIP) in Sec. IV-B that is compatible 
with MIP solvers such as CPFEX [20]. Finally, it is shown in 
Sec. IV-C that the relationship Ik = fk{Pk) is linear. 

A- The Visit Lifetime Maximization Problem 

In this section we propose a continuous, non-linear opti¬ 
mization problem for maximizing the visit lifetime Ik when 
given a visit expenditure pk, assuming the network model 
in Sec. III-A. By definition the visit number fc > 1. We 
let a = [oi 02 ... where at is the energy allocated 

to each location i € U- We express our budget Bk{a) in 
terms of a. As we will show below, the energy vector can 
be used to determine the node hardware expenditure, Xk, 
and the energy expenditure, Yk- The sum of these costs with 
the labor expenditure, Zk, is equal to the overall budget. 
By constraining the budget Bk{a) to the visit expenditure 
Pk, we can optimize the node hardware, energy, and labor 
expenditures without causing the budget to exceed pk- We also 
provide the power P (i) in Watts consumed at each location i 
in terms of the information flow rate , the rate at which each 
location i sends data to other locations j- After discussing the 
expressions for Bk{a) and P (i), we propose the optimization 
problem itself. 

The budget is a function of the location energy allocation 
vector a. The energy expenditure Yk = which 

requires the energy at in Joules allocated to each location 
i G U, as well as the cost a of a single Joule. To determine 
whether or not a node is required at location i, we notice that 
location i requires a node if O-i 1> 0. We use an indicator 
function to tell us whether location i needs a node based 
on the value of a^. To obtain a continuous formulation, the 
exponential function (1 — exp(—Coi)) can be used as an 
indicator function. The constant C is a large number such that 
the function is approximately 1 when ai > 0, and 0 otherwise. 

We assume that all nodes have identical hardware, and are 
purchased at a price of /3 financial units per node during the 
initial visit. We therefore let Xk = P (1 — exp(—Coi)) 
when k = 1, and Xk = 0 when fc > 2. In this case, all nodes 
are assumed to have the same cost (3 and the same rate of 
failure A, regardless of whether they are sensor, relay, or sink 
nodes. 

As we assume that node hardware is identical and that 1- 
connectivity is sufficient, each node has the same failure rate 
A, and a single node failure potentially results in network 
failure. The overall Mean Time Between Failures (MTBF) of 
the network is therefore 


Him ^ (1 - exp(-Cai)) AXi' 
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By optimizing the node hardware expenditure Xi during the 
first visit, we are also maximizing the MTBF w. We do not 
consider placing redundant nodes in this section, as that is 
covered in the method for optimizing node hardware costs and 
the approximate NPC of unscheduled payments in Sec. V-D. 

We assume that the labor expenditure Zk during visit k is 
known in advance. Later in this section we show that, because 
we optimally allocate energy to the nodes, every location 
with a node will be visited when maintenance is performed. 
This means we know which locations to visit in advance, and 
can estimate the time required to travel between the different 
locations, and therefore the time spent and cost of performing 
labor. When factoring in node hardware, energy, and labor 
costs, our budget can be expressed as 

B (a) = I “ exp(-Cai)) + aa^) + Zk k = l 

^ + > 2 . 

(5) 

The power P (i) in Watts consumed by node i gU can be 
expressed in terms of the information flow rates rij. Location 
i sends data to location j at a rate of rij bits/s. We let the 
sets Mi, Ti, and Ui represent the set of nodes in M, T, and U 
respectively that are within communication range of node i. 
The energy consumed by node i when transmitting to another 
node j is Ti {dij) J/bit, where dij refers to the distance between 
the nodes in m. Node i consumes Qi J/bit while receiving data. 
Sensor i consumes Si J/bit collecting data at a rate of gi bits/s. 
The expression for the power in Watts consumed by a node is 

-P (0 = X! + X] 

j€Ui j&XiUTi 

We can use (5) and (6) to formulate the lifetime maxi¬ 
mization problem as a non-linear, non-convex, yet continuous 
optimization problem 


max Ik 

a.i,rij,lk G R+ 



(7a) 

S.t. ^ Tp + g. 

= X M 

Vi e A/'ur 

(7b) 

j&AfiUTi 

hP (i) = ai 

j&Ui 

\/iGU 

(7c) 

Bk{a) = Pk- 



(7d) 


The energy allocated to node i, the information flow rate 
Vij between node i and node j, and the visit lifetime Ik are all 
optimization variables that exist in the non-negative reals. The 
elements of vector a are the energy allocation values a^. Flow 
constraint (7b) ensures each node does not transmit more data 
per second than it receives from others or collects through 
sensing. The sensor data generation rate gi > 0,Vi G T 
because sensors generate data, while gi = QMi G M because 
relays do not. Energy constraint (7c) ensures each node i uses 
all of the ai Joules allocated to it; we explain later in this 
section why this is an equality constraint. Budget constraint 
(7d) lets us find the optimal node hardware, energy, and labor 
expenditures resulting from adjustments to a, while ensuring 
that the budget does not exceed the visit expenditure pk- 
The left-hand sides of (7c) and (7d) are non-linear; as both 
constraints are equality constraints, (7) is non-convex. 


The energy constraint (7c) is typically defined as an in¬ 
equality constraint (for example, in [3], [4], and [21]). To 
demonstrate why representing (7c) as an equality constraint 
is valid when maximizing lifetime and optimally allocating 
energy, we will prove by contradiction that at optimality the 
left-hand side of (7c) must be equivalent to the right-hand 
side. Suppose that at the optimal lifetime l^P (i) < ai for 
a node i. This implies that 

1) At optimal lifetime one or more nodes have exhausted 
their supply of energy. 

2) Node i has (a^ — l^P (/)) J of energy that it has not yet 
spent. 

If, when adding or replacing batteries, a portion of the spare 
energy in node i had instead been allocated to the nodes 
with no energy, a lifetime longer than would have been 
achieved. This means that is not in fact optimal if 3i G 
U : IIP {i) < Qi- Therefore is only optimal as long as 
l^P (i) = UiMi G U. This further implies that all locations 
with nodes will be visited when replacing batteries, as all 
nodes will have exhausted their energy supplies at the optimal 
lifetime 

The same argument can be used to justify making the budget 
constraint (7d) an equality constraint. We can prove by contra¬ 
diction that, at the optimal lifetime the budget i?fe(a) = pk- 
Suppose that at the optimal lifetime Bk{a) < Pk, meaning 
that we have extra money {pk — Bk{a)) that has not been 
spent. Additional energy could have been purchased with this 
money, increasing ai,Vi G U until Bk{a) = pk- From (7b) 
we know that the flow rates will not increase with energy 
ap. the rate gi at which data is generated by sensor i is fixed, 
therefore according to the energy constraint (7c) lifetime 
must increase. This implies that is not in fact optimal, 
meaning an optimal lifetime Z^ requires that Bk{a) = pk- 

B. Formulation of the Mixed Integer Program 

The optimization problem given in (7) is non-convex, mak¬ 
ing it difficult for a solver to find a globally optimal solution. 
Even though (7) is non-convex, it can be re-written as a Mixed 
Integer Program (MIP) and can therefore be solved with the 
robust branch-and-bound algorithms and heuristics available 
in commercial MIP solvers. The constraints will all be made 
linear, making the problem compatible with MIP solvers such 
as CPLEX [20]. 

Multiplying (7b) by Ik allows us to express the flow of 
data between node i and node j as bij bits instead of the rate 
Tij bits/s. The left-hand side of (7c) is re-written in terms of 
bij. To turn the continuous formulation (7) into a MIP, we 
replace the indicator function with the binary optimization 
variable Xi that is 1 when the node either exists at or will 
be added to candidate relay location i, and 0 otherwise. With 
these changes in mind, the optimization problem (7) can be 
rewritten as 

max Ik (8a) 

ai,bij ,1}^ G 

s.t. ^ bji + gih = X Vf e A/’U T 

jGMiDTi jeUi 


(8b) 
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^ ^ “t" ^ ^ Qibji 

3&ii jeJ^iUTi 

+ SiQilk = 0 ,i \/i ^lA 

(8c) 

{Pxi (1 - i/[fc - 2]) + aai) 

i&A 

+ Zk=Pk (8d) 

Qi < Dxi 'ii ^lA 

(8e) 

Xi = l Vi £ TU <S, 

(8f) 

where D is a constant greater than or equal to the largest 
possible battery capacity. The flow, energy, and budget con¬ 
straints (8b), (8c), and (8d) serve the same purposes as 
their counterparts (7b), (7c), and (7d) respectively. Note that 
H\k — 2] represents the Heaviside step function, making 
I3xi (1 — i7[fc — 2]) = j3xi when fc = 1, and 0 otherwise. 
This is done because, as discussed in Sec. IV-A, nodes are 
only purchased during the initial visit. We ensure that energy 
is only added to location i if it has a node via (8e). Constraint 
(8f) ensures that sensor and sink nodes are always placed on 
the field. The objective and all the constraints in (8) are in a 
linear form; note that the Heaviside function is solved prior 
to optimization, meaning that (8d) is linear. The optimization 
variables are either real numbers or integers, making (8) 
compatible with commercial MIP optimizers. 


V. Framework for Optimizing Net Present Cost 

When Ik is maximized for a given pk, we established in 
Sec. IV that the relationship Ik = fk{pk) is linear. In this 
section we utilize this relationship to minimize the overall 
WSN NPC expressed in (3). 

Sec. V-A provides a general non-linear, non-convex NPC 
minimization formulation that is applicable to any Ik = fk{Pk) 
lifetime function, linear or not. This generalized formulation 
is non-convex and difficult to solve optimally. However, in 
Sec. V-B we demonstrate that when the linear lifetime function 
h = fk{Pk) defined in Sec. IV is used, NPC minimization 
for a fixed number of visits K is convex. Assuming that a 
maximum of K visits may be made, we provide a 0{K^) 
algorithm for minimizing the NPC with a linear lifetime 
function. In Sec. V-C we show that equally spacing the visits 
can be a good rule of thumb when fk{Pk) is linear. 

The NPC minimization framework assumes that the net¬ 
work’s Mean Time Between Failures (MTBF) is known and 
constant. The network’s MTBF and therefore the NPC of 
unscheduled payments might depend, however, on the cost 
of node hardware: using more expensive nodes or adding re¬ 
dundant nodes in stand-by to each node location may improve 
the MTBF. A method of balancing the initial cost of node 
hardware with the approximate NPC of unscheduled payments 
is provided in Sec. V-D. 

A. General NPC Minimization Framework 


C. Lifetime Function Derivation 

To minimize the overall NPC in Sec. V, it is necessary to 
express the results of the Ik maximization in this section in 
terms of the Ik = fk{Pk) function in (2). By doing so, we 
demonstrate that fk{Pk) is a linear function. 

After running the optimization described in the previous 
sections, the maximum visit lifetime and the optimal energy 
a* at each node i that achieves it, have been found. The budget 
function (5) can therefore be rewritten as 


Bkia*)=Xk+allp + Zk, (9) 


where a* is a vector whose element i is a*, and p — 
Thi&A ^ (*) represents the overall network power consumption 
after the optimal information flow rates r*j have been deter¬ 
mined. Constraint (7c) allows us to substitute a* with I'^P (i), 
therefore allp = (*)) = The node 

hardware expenditure Xk = ~ e^P(~C'a*)) when 

fc = 1, and Xk = 0 otherwise. 

Rearranging (9) and letting Bk{a*) = pk as in (7d) gives 


l*k = fkiPk) 



Pk > Xk + Zk 
Pk < Xk + Zk- 


( 10 ) 


The terms Xk, Zk, a, and p do not change with pk- From (7b), 
the flow rates depend on the sensor data generation rate pi, 
which is constant. From (6) we can see that P (i) and therefore 
p only change with and pi, so p is not a function of pk- 
The flow rates and therefore the locations requiring nodes do 
not change with pk, so Xk does not change with pk- Both Zk 
and a are constants. As a result, increases linearly with pk- 


Our goal is to minimize the NPC in (3) while ensuring 
that the network remains operational for exactly L years. 
As discussed in Sec. III-B, the relationship between visit 
expenditure pk and visit lifetime Ik of visit fc is the function 
fk{Pk)- By finding the optimum visit expenditures we And 
the optimum visit lifetimes = fk{Pk) ^nd thus the optimum 
spacing between visits. 

To minimize the NPC, we require upper bounds on the visit 
expenditures and the number of visits. The maximum number 
of visits is K, meaning K < K. For example, if the network 
operator could not feasibly visit the network more than once 
per month, then K = 12L. The maximum visit expenditure 
for visit fc is denoted pk, meaning pk < Pk,'^k £ [1, K], 

The NPC minimization problem for when K > 2 can be 
written as 


K 


P1 + J2 


Pk 


mm x-j- ’ / k-L 

pi...p*eR+ ^ (1 

xfc,xfc+ie{o,i} 

LL/C.J 


E 


^ (i+t;)" 

s-i- ^fk{Pk)=L 


K 


k=l 

Pk < XkPk 


Xk+1 Xk 


(11a) 


(lib) 

vfc £ [i,k] 

(lie) 

Vfc £ [i,k-i] 

(lid) 
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fk{Pk) > 0 


VfcG 

(lie) 


where the binary value Xk represents whether or not a network 
operator performs visit k, and w > 0 represents the rate at 
which money earns interest. The constants uj and cr represent 
the network’s MTBF and the cost of network failure respec¬ 
tively. Note that when K = 1, finding the NPC is trivial: 
the visit expenditure pi that achieves the operational lifetime 
L is chosen. The objective function (11a) is the NPC, while 
the lifetime constraint (lib) ensures that the visit expenditures 
provide the required operational lifetime of L. Constraint (11c) 
makes Pk = 0 when visit k is not made, while guaranteeing 
that Pk does not exceed the maximum expenditure amount 
Pk when a visit is made. Constraint (lid) is used to find 
the optimum number of visits: if Xi = 0 for visit i, then 
Pj — G [i,K], making the optimum number 

of visits K* = J2k=i^k- Constraint (lie) ensures that visit 
lifetime can never be negative. 

The NPC minimization problem (11) is a non-linear mixed 
integer program. If fk{pk) is non-linear, then (11) is non- 
convex. Suppose that iT > 2, q = [pi p 2 ■ • ■ Pj^], and r(q) 
is (11a). It can be shown that 


c)^r(q) £)^r(q) 

dptdpp^ dpj^dpr 


In(l-r^)) dfijpi) 

0 


i ^ k 
i = k. 


We can see that the element at row K and column K of 
the Hessian of r(q) is zero, and the other elements in row 
k and in column k are non-zero. The Hessian is therefore 
not positive semi-definite, meaning that the objective function 
(11a) is non-convex. When iC = 1, finding the NPC is trivial 
because pi is set to the value that ensures /i(pi) equals 
operational lifetime L. 


B. NPC Minimization with a Linear Lifetime Function 

Based on our proof of linearity of the lifetime function 
fk{Pk) in Sec. IV, this section will demonstrate that the NPC 
minimization problem is convex for a fixed number of visits 
K. We provide an algorithm to determine the optimal number 
of visits in 0{k^) time, for the maximal number of visits k. 

As shown in Sec. IV-C, we can assume the lifetime function 
has the form 


optimal number of visits K* when given the maximum number 
of visits possible K. Doing so allows a network operator 
to find the optimum visit expenditures, visit lifetimes, and 
number of visits. 

We can derive an equation for NPC when K visits are made 
and fki^k) = knpk + bk when pk > —bk/rri. Using the fact 
that J2k=i fk{Pk) = L, and letting p = [pi p 2 • ■ • Pk], we 
can re-write the NPC (3) when K > 2 as 


9{P) 


L-bi- Y.k=2 + hfc) 

m 

K 

_ Pk _ 

^ (1 -f 

IL/uij 

y _^_ 

r> = 1 ' ' 


Note that the first term represents the initial visit expenditure 
Pi. When K = 1, the initial expenditure pi is set to the value 
that ensures an operational lifetime L. 

We can find the visit expenditure vector p = [p 2 ■ • ■ Pk] 
that minimizes the NPC ^(p) for K >2 using the Lagrangian 


K 

u) = 5 (p) - y Uk{pk + bk/m), 

k=2 


where u = [u 2 ... uk]- When solving a minimization prob¬ 
lem with inequality constraints, the slack variables Uk are op¬ 
timized to satisfy the Karush-Kuhn-Tucker (KKT) conditions. 
The complementary slackness KKT condition implies that 
Pk = —bk/m when Uk > 0, and Uk = 0 when pk > —bk/m. 

We can use the gradient of the Lagrangian to derive an 
equation for finding pk+i from pk- The stationarity KKT 
condition requires V£(p,u) = 0, where VC(p,u) = 

r 3£(p,u) a£(p,u) a£(p,u) 9t:(p,u) i c„Ufraprin(i 

i Opi dpK dui duK V Subtracting 

_ 9£(p,u) k < K — 1, yields 

apk+i apk ■> — ’ 


Pfc +1 = q{pk) 




TO ln(l-f r;)(l-b mpfc+i &fc+i ’ 


where 


9{Pk) = 


(l + v) 


mpk+bk 


- 1 


TO In (1 -b v) 

The visit expenditure Pk+i can be found using 


fk{Pk) 


mpk + bk pk>-bk/m 
0 pk< -bk/m, 


Pk+i 


q{Pk) q{Pk) > -bk+i/m 
-bk+i/m q{pk) < -bk+i/m. 


( 12 ) 


where to > 0. Note that the slope to is not a function of k. A 
constant slope to for all visits implies that the cost of energy, 
and the rate that the network consumes energy, do not change 
over time. Using (10) from Sec. IV-C, the slope m = -^ 
and y-intercept bk = — ■ As with the general NPC 

minimization problem, we find the optimum visit lifetimes 
Ik ~ fk{Pk) finding the optimum visit expenditures 

Pi- 

In this section we first provide a method for finding the visit 
expenditures pk for fc G [l,Ar] to minimize NPC when K is 
fixed. We then provide a 0{k^) algorithm for determining the 


Keeping in mind that —bk+i/m > 0, if q{pk) > —bk+i/m, 
then Pk > -bk/m so Uk = 0. As pk+i > q{Pk) > -bk+i/m, 
Uk+i = 0 so pk+i = q{pk)- When q{pk) < -bk+i/m, 
then Uk+i > Uk to ensure that Pk+i > —bk+i/m-, to 
satisfy the complementary slackness KKT condition, making 
Uk+i > Uk>0 forces pk+i = -bk+i/m. 

The KKT conditions are maintained by (12), meaning it can 
be used to find a local optimum solution. By ensuring that 
Uk = 0 when pk > —bk/m, and pk = —bk/m when Uk > 0, 
it satisfies the complimentary slackness condition. The primal 
and dual feasibility conditions are satisfied by ensuring that 
















Pk > —bk/rn and Uk > 0 respectively. To derive (12), the 
stationarity condition V£(p,u) = 0 was assumed. 

When q{pk) < —bk+i/m and Pk+i = —bk+i/m, note that 
no energy is being added to the network. Visit k + 1, and every 
subsequent visit, will have a visit lifetime of 0. If this occurs, 
the chosen K value is too high and is not the optimal K* 
value. 

Before discussing our method of using (12) to minimize 
g(p), we first prove that the local minimum of g{p) is also 
the global minimum by showing that ^(p) is convex when 
K >2. We will do so by letting 

T /T \ _ _Pfe_ 


and proving that h{k) is a convex function with respect to the 
visit expenditures pk- If h{k) is convex, then ^(p) is the sum 
of convex functions and is therefore itself convex. The second 
partial derivative of h{k) is 


d^hjk) 

dpidpj 


' 2m ln(l+'u)+m^ ln^(l+i;)pfe 

m ln(l+i?)+m^ ln^(l+'i;)pfc 
< (l+„)r'-E^=fc(™Pn+l.n) 

ln^(l+«)pfc 

(l+^)r'-E;)=fc('pp7.+i>n) 

.0 


i = k A j = k 

{i = k Aj > k) 

V ( j = k Ai > k) 
i > k A j > k 

i < ky j < k. 


By definition pk > —bk/m > 0, ensuring that the elements 
of the Hessian H of h{k) are all non-negative, and that the 
elements of the visit expenditure vector p are non-negative 
real numbers. This means pHpT > 0, so H is a positive 
semi-definite matrix. As H is positive semi-definite, h{k) is 
a convex function, so g{p) is the sum of convex functions 
and is therefore convex itself. The local optimum found when 
minimizing NPC for K visits is consequently the global 
optimum. 

Alg. 1 can be applied to find the global minimum of 
p(p). It determines the optimal value of p 2 , then applies 
(12) to find the remaining visit payments. The function 
FlNDR00T(/(a;)) finds a; > 0 such that f{x) = 0. The 
derivative of the NPC’s Lagrangian £(p,u) is ex¬ 

pressed in terms of p 2 and K by dLagrangian(p 2 , K). As 
the local minimum of g(p) is its global minimum, calling 

FindRoot(dLagrangian(p 2 , K)) finds the value of p 2 that 
minimizes p(p). The function NEXTP(pfc, AT) uses p 2 and (12) 
to calculate the visit payments p^,p 4 ,... ,px- 


Algorithm 1 Optimizing the visit expenditures for a fixed 
number of visits K. _ 

// Returns the visit expenditure vector p that achieves the 
// optimal NPC if Ff > 2 visits are performed, 
function OPTIMALPAYMENTS(Ar) 

// Find the optimal value of p 2 

P2 A- FlNDROOT(DLAGRANGIAN(p2,f^)) 

// Find the remaining payment values 
for /c ^ 2 to {K — 1) do 
Pk+l ^ NEXTP(pk,K) 

_ ^ L-bi-J2k=2('^Pk+bk) 

Pi ^- m - 

return p 


We use Alg. 2 to determine the number of visits K that 
optimize the NPC. It finds the NPC for every K G [1, FT] using 
Alg. 1, and chooses the values K* and p* that minimize the 
NPC. The function NPC(p, K) uses (3) to calculate the NPC. 
Alg. 2 is bounded by 0{K^) operations. In order to express 
Pk in terms of p 2 , (12) must be applied k — 2 times. The 
number of operations in the equation for after all pk 

values are expressed in terms of p 2 , are bounded by 0{K^). 
Expressing K—2 values in terms of p 2 requires applying (12) 
a total of {K — 1){K — 2)/2 times, leading to the 0{K‘^) 
bound. Once p 2 is found, the number of operations to find 
the remaining K — 1 visit expenditure values are bounded by 
0{K). As K G [1, K], we have to repeat these steps K times, 
meaning that finding the optimal number of visits K* and the 
corresponding optimal NPC g(p*) is an 0{K^) algorithm. 


Algorithm 2 Finding the optimal number of visits and NPC 
value for up to K visits. 


function OPTIMALNUMBEROFVlSITS(iT) 


K* G- 


L-bi 

rr 

1 


// Try every value of K 
for AT ^ 2 to A' do 

p ^ OPTIMALPAYMENTS(Ar) 
g G- NPC(p,A:) 

if g < g* then 

9* ^ 9 
K* ^ K 

return AT* 


C. Equal Visit Lifetime Approximation 

Equally spacing the visits apart can be a good rule of thumb 
when the lifetime function fk{pk) is linear as in Sec. V-B, 
and \bk\ << l,Vfc G [2, A']. In the previous sections the visit 
lifetimes could be non-uniform; they could vary for each visit 
k depending on the visit expenditure pk- In this section we 
prove that visit lifetimes Ik are approximately equal as long 
as V « 1 and |6fc| << l,Vfc € [2, AT]. 

We assume that the value of K is chosen such that, Vfc G 
[1, AT], q{pk) > -bk+ilm, so pk+i = q{pk)- We can use (12) 
to express the lifetime Ik+i of visit /c -f 1 in terms of the visit 
expenditure pk of visit k through 

+ _i 

IkGl — -;— 1-f Ofc-l-1- 

In (1 -f v) 

We can write = exp[ln(l-|-u)(mpfc-|-6fe)], and as 

long as ln(l-|-u)(mpfe-|-6fc) << 1, exp[ln(l-|-u)(mpfe-|-hfc)] 

1 + ln(l + v){mpk + bk), so after simplification 

^fc-i-i ~ GbiPk + bk + bk+i, 

which can be simplified to 

^fc-i-i ~ (13) 

As long as v « 1 and \bk\ « l,Vfc G [2, AT], evenly 
spacing the lifetimes is a near-optimal approach. If u << 1 
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and \bk\ << 1, then it is likely that lii{l+v){mpk + bk) << 1, 
which is required by the approximation. The magnitude of bi 
does not have to be small; in fact, chances are it will not be 
small due to the node hardware costs of the initial visit. We 
call this the Equal Visit Lifetime (EVL) approximation. 


D. Minimizing the Initial Cost of Node Hardware and the 
Approximate NEC of Unscheduled Maintenance 

To reduce the NPC of unscheduled payments, a network 
operator could purchase more robust nodes or add redundant 
nodes in stand-by to node locations in order to increase the 
network’s reliability. This, however, increases the initial cost of 
node hardware. By increasing the initial cost of node hardware 
Xi, the network operator can increase the spacing between the 
failure times ti.. .tp of the network and reduce the number of 
failures F that occur over the network’s operational lifetime. 
The NPC (1) will decrease if the reduction in ^^=1 (i+vf'^ 

exceeds the increase in Xi. When Xi and X]n=i (i+vy-n 
have been found, the NPC can be minimized using (11) if the 
lifetime function is non-linear, or using the methods provided 
in Sec. V-B and Sec. V-C if the lifetime function is linear. 
This section provides a method for balancing Xi, tn, and F 
to minimize the cost of node hardware and the approximate 
NPC of unscheduled payments in (3), where u: is the network’s 
MTBF, tn = nuj, and F = [T/wJ. 

We assume that the network operator is given a set of 
choices %, with each choice i costing jdi financial units per 
node location while having a failure rate of Xi per node 
location. Choices could represent different node hardware, 
each with their own cost and failure rate. Alternatively, the 
choices could also include different levels of redundancy. 
Assuming Poisson failures, putting G nodes at a location 
would reduce the failure rate by a factor of G while increasing 
the cost by a factor of G. 

We assume that node hardware is only added during the 
initial deployment, that all nodes at a single node location 
failing leads to a network failure, that the network operator 
knows the number of nodes N that are required by the 
network, that every node has the same type of node hardware, 
and that every node location is given the same number of 
redundant nodes in stand-by. The approximate cost of each 
choice with respect to the initial payment and unscheduled 
payments can then be written as 


[LN\i\ 

N^3^+ 

n—1 


O' + Pi 
(1 ’ 


(14) 


where cr is the cost of performing the repair. To find the 
optimum choice i*, the operator can try every choice and 
pick the one that minimizes (14). The number of operations 
required are bounded by LNXf). 


VI. Numerical Results 

In this section we provide numerical results using the 
frameworks from Sec. IV and Sec. V. Sec. VTA presents 
node and network assumptions based on a gas monitoring 
scenario. We analyse the effectiveness of our method for 



Fig. 2. A possible area for monitoring in our scenario. Map data and image 
©2014 Google. 


balancing initial node hardware costs with the approximate 
NPC of unscheduled payments in Sec. VTB. In Sec. VTC 
we demonstrate that NPC minimization is effective at reduc¬ 
ing costs and that equally spacing visits is a near-optimal 
strategy. Sec. VTD demonstrates that NPC minimization is 
most effective when node hardware and labor expenditures 
are low and a large portion of the budget is spent on energy. 
Sec. VTE examines how adjusting the number of visits K, 
the operational lifetime L, the interest rate v, the network’s 
MTBF w, and other parameters related to the cost of node 
hardware and labor affects the portion of the NPC allocated 
to node hardware, failure repair, energy, and labor. 

The percent savings is used in this section to measure 
the effectiveness of NPC minimization when compared to a 
network design based on a single payment at the start of 
the network lifetime. The minimized NPC value is the cost 
g{[pi ■ ■ -Pk]) of running a network for L years by performing 
NPC minimization over K visits. The single payment bench¬ 
mark is the cost, s, of providing each node with enough energy 
to last L years at the time of initial deployment; note that s is 
the NPC when K = \. The percent savings is defined as the 
percent difference between these two values such that 

\s - g{[pi .. .pk])\ 
s 

In this section we use the dollar symbol $ to represent a 
financial unit. The NPC minimization framework does not 
assume a specific currency; any currency can been used as 
long as it is consistent for all costs. Prices given for node 
hardware, energy, and labor reflect their value in United States 
dollars at the time this paper was written. 

A. Node and Network Assumptions 

We consider applying the network and cost models from 
Sec. Ill to a gas monitoring scenario, where a company could 
install a WSN to continuously monitor the concentration of 
hazardous gas in the air at a storage site. The satellite image 
in Fig. 2 provides an aerial view of one area where such a 
network could be built. When finding the parameters for our 
numerical results, we assume that sensor nodes are the circles 
positioned on top of the tanks in the picture. The sink node is 
located in the middle of the square. To connect disjoint sensors 
to the sink, we populate the candidate relay location set, M, 
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with locations spaced along a minimum Steiner tree generated 
by the GeoSteiner [22] tool using the process described in [23]. 

In our scenario, we assume that the relay and sensor 
nodes are similar to Iris motes [2] and that TGS 825 gas 
sensors [24] are used. The gas sensor is heated for about 60 s 
prior to each measurement, meaning 39.6 J are consumed per 
measurement; the measurements occur every 5 minutes. The 
energy consumed by the RF230 transceiver and ATmegal281 
microcontroller present on Iris motes while receiving and 
transmitting data are taken from their datasheets [25] and 
[26] respectively. We also assume they are connected to a 3 V 
source, that the RP230’s data rate is 250kbit/s, and that nodes 
transmit data in 64bit packets. Note that in general, sensor 
nodes measuring gas concentration require a large amount of 
energy due to the heating elements in the gas sensors, and 
consume significantly more energy than the relay nodes in the 
network. The results in this section, however, do not rely on 
sensor nodes consuming more energy than relay nodes. Similar 
results would be found in a scenario where relays consume 
more energy than sensors, such as a video surveillance network 
where relays may have to forward large amounts of data 
frequently, as long as the network consumes power at a similar 
rate. 

We minimize the power Ti{dij) consumed by the RF230 on 
node i by adjusting its transmit power when sending data to a 
node j that is dij m away. The transmit signal power P™ W 
required at node i to ensure a receive signal power W at 
a destination node j is calculated using the Friis equation 


pTX ^ 

ij 


P^^{AT:dof fd., 
G,G,A2 


(15) 


where dg = 1 m is the distance from the antenna to the edge of 
the near field, Gi = G^ = 1.5 is the antenna gain of a dipole 
antenna, A = 125 mm is the wavelength of the signal, and 7 = 
4 is the path loss exponent. We set the receive signal power 
P^ = -101 dBm, which is the receiver sensitivity of the 
Iris mote. We calculate the 16 discrete RF230 transmit signal 
powers with a quadratic interpolation of the transmit current 
consumption values for different transmit signal powers given 
in its datasheet, which are then used with (15) to determine 
the power consumed by a node i when transmitting to another 
node j that is dij m away. 

As stated in Sec. IV-A, we assume that the node hardware 
cost consists only of purchasing the nodes for initial deploy¬ 
ment. We adjust the cost of a single node, /?, between $10 and 
$ 100 ; such a range covers cases where an operator builds 
nodes independently and where an operator purchases them 
from a manufacturer. By running ( 8 ) on the network in Fig. 2, 
we find that about 150 nodes in total must be purchased. 

In Sec. IV-A we discuss that the portion of an expenditure 
dedicated to energy depends on the cost of each Joule and 
the amount of Joules required by the network. The cost of a 
Joule a when using alkaline D-Cell batteries is approximately 
20 p$/J, but is 66 p$/J for Lithium D-Cell batteries, which 
are lighter, have a higher capacity, and are better suited 
for extreme climates. The network described in Sec. VI-A 
consumes energy at a rate p of approximately 6.2 W while 
running. 


We assume that the price of labor Zk for each visit k 
does not change over time; that is Vfc G [1,A"] ■ = C 

in dollars. To find the labor cost we ran the Traveling 
Salesman Problem [27] on the network in Fig. 2; when 
assuming a walking speed of about 1.4m/s [28], it takes 
approximately 2 hours for someone to visit all of the nodes. 
With a wage of 20$/h, and assuming 60s per node is taken 
for swapping batteries, we assume ( is minimum $140. The 
labor cost C, would likely be higher for networks with nodes 
that are difficult to reach, networks that are in remote areas, 
or networks in areas that are difficult to traverse. We assume 
that a visit could cost up to $1000 in such scenarios, due to 
higher wages, extra time spent at the site, and transportation 
costs in order to reach the site. 

We assume that all nodes have identical failure rates. While 
we were unable to find failure rate data on commonly used 
motes, such as the Iris mote, we were able to obtain wired gas 
monitor failure rates from Draeger [29]. The DraegerSense IR 
has one of the best failure rates of about 0.5 pfailures/h. While 
Draeger sensors are not wireless, the CC2420 transceiver has 
a failure rate in the order of 1.9 nfailures/h [30]. We assume 
that a high-quality node with a gas sensor could have a slightly 
higher failure rate of 0.75 pfailures/h. We assume lower-quality 
nodes have failure rates of 10 pfailures/h, higher than that of 
Draeger Polytron 7000 gas monitors with electrochemical sen¬ 
sors. The high-quality and low-quality node hardware failure 
rates correspond to network MTBF ui of about 1 year and 
1 month respectively. 

We set the cost of each repair to $1000. While it takes less 
time to replace a single sensor than to restore energy to all 
N = 150 nodes, we cannot predict when failures occur, so 
unlike replenishing energy, there will be downtime between 
the beginning of the failure and its repair. The repair cost value 
accounts for the network being unable to perform its function 
during this downtime, as well as the cost of having some type 
of emergency responder available to repair the network. 

We placed limits on the operational lifetime of the network 
and the interest rate when minimizing NPC. Operational 
lifetime, L, ranges from 1 to 10 years. Based on an average 
stock market rate of return of 0.083 over the last 114 years in 
developed countries [31], the interest rate is assumed to vary 
between 0.01 and 0.1 to capture below and above average rates 
of return. 

In the following sections, results are generated by assuming 
a default set of parameters. We assume the cost /3 of each node 
is $10, 150 nodes are present in the network, the cost a of a 
single Joule is 20 p$/J, the network consumes energy at a rate 
p of 6.2 W, the cost of labor is 1000$/visit, the MTBF ui 
is 1 year, and the cost a of each failure repair is $1000. By 
default, the operational lifetime L is 10 years, and the interest 
rate u is 0.1. All results in the following sections are generated 
using these default values unless stated otherwise. 


B. Analysis of Minimizing the Initial Cost of Node Hardware 
and the Approximate NPC of Unscheduled Maintenance 

In this section we demonstrate how the initial cost of node 
hardware and the approximate NPC of unscheduled payments 
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Number of Nodes per Node Location 

Fig. 3. Normalized value of (14) for various numbers of nodes placed at 
each node location in the network. 

can be minimized using the method described in Sec. V-D. 
We assume that we have cheap, low quality nodes, and we 
want to find the optimal number of redundant nodes to place 
in stand-by at every node location in the network, assuming 
that the network fails if all the nodes at any one node 
location fail. We demonstrate that the method in Sec. V-D can 
significantly reduce the approximate NPC of initial hardware 
and unscheduled payments. 

Our choice set % consists of placing 1 to 10 nodes at each 
node location, where one node per location is active and the 
rest are in stand-by. Each node costs $10 and has a failure rate 
of lOpfailures/h. We assume Poisson failures, so putting G 
nodes at a location reduces the failure rate by a factor of G 
and increases the cost by a factor of G. 

Fig. 3 shows the value of (14) for each number of nodes 
per node location tested, normalized to only deploying one 
node at each location. The lower bar represents the value of 
node hardware, while the upper bar is the approximate NPC 
of unscheduled payments. 

We can see in Fig. 3 that the method in Sec. V-D can 
significantly reduce the approximate NPC of node hardware 
and unscheduled payments. Node hardware costs increase 
linearly, but for small numbers of redundant nodes the ap¬ 
proximate NPC of unscheduled payments decreases rapidly. 
After a certain point, for example in Fig. 3 after 7 nodes have 
been placed at each node location, savings from additional 
redundancy no longer outpace additional hardware costs. 

C. Comparison of Optimization Strategies 

In this section we compare the Equal Visit Fifetime (EVE) 
approximation to finding the Optimal NPC (ONPC). The 
lifetime function is linear, so the ONPC is found using the 
approach described in Sec. V-B. Fig. 4 shows how the NPC 
found with both methods depends on the total number of visits 
K and the labor cost ( in $/visit. NPC is normalized to the 
cost of a network with an operational lifetime of L = 10 years 
when K = I and C =1000 S/visit. The value of ^ written 
underneath each pair of EVE and ONPC lines was used to 
generate that pair of lines. The minimum ONPC value over 
all K for a given value of ( is indicated by the squares. 



Fig. 4. NPC for various numbers of visits when the cost of labor vaiies 
between 140$/visit and 1000$/visit. 

The squares in Fig. 4 show that as the labor cost ( decreases, 
the minimum NPC over all numbers of visits K decreases, 
while the optimal number of visits required to achieve that 
NPC increases. Decreasing the labor cost means that more 
visits can be performed before the cost of performing labor 
exceeds the savings of NPC minimization. More visits provide 
additional opportunities to exploit interest rates to purchase 
cheaper energy, reducing the money spent on energy to achieve 
a given operational lifetime L, lowering the overall NPC of 
the network. 

The rate of change of the EVE and ONPC lines are greater 
for higher values of C,. The higher the labor cost, the greater 
the total cost of labor for a given number of visits K, and 
the greater the rate of change of the lines generated with 
that labor cost. When K < 5 visits, NPC decreases rapidly 
with increasing K because additional visits provide more 
opportunities to use interest rates to reduce the cost of energy. 
Eventually increasing the number of visits becomes ineffective 
at reducing energy costs. For instance, the slopes of the 
EVE and ONPC lines are positive when K > 5 visits and 

C, =1000$/visit, and the slopes of the EVE and ONPC lines 
are flat when K > 10 visits and ( =140$/visit. 

The optimal NPCs calculated by the EVE approximation are 
within 0.01% of the NPCs calculated using ONPC. For large 
values of K, however, the NPC found by EVE deviates from 
that of ONPC when C is high, as shown in Fig. 4. In Sec. V-C 
we assume that |6fc| << 1; however, the |6fc| values are not 
<< 1 for high (. This can cause the difference between the 
Pk values calculated by EVE and ONPC to be relatively large, 
noticeably affecting the NPC when K is also large. When the 
cost of labor is high, the optimal number of visits tends to be 
low, meaning that even though the error is present, it is not 
significant in the range of K where NPC is minimized. 

D. Analysis of NPC Minimization Parameters 

In this section we examine how the performance of NPC 
minimization is affected by changes in the parameters de¬ 
scribed in Sec. VTA. To determine the situations where NPC 
minimization is most effective under our network model, we 
will find the parameters that have the biggest impact on the 
optimal NPC. The parameters considered are the cost [3 of 
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Fig. 5. Percent savings of NPC minimization compared to performing a 
single visit, for various energy payment rates 0 at different node hardware /3 
costs, different labor costs, and different MTBFs uj. 

node hardware, the cost a of a Joule, the rate p of energy 
consumption, the cost of labor, the interest rate v, the 
network’s operational lifetime L, and its MTBF uj. 

Two parameters that influence the return earned on an 
investment are L and v, NPC minimization uses the returns 
provided by interest rates to reduce future energy expenditures. 
If either L or u is low, the percent savings is less than 5%, 
meaning that NPC minimization is ineffective. When L is low, 
little time is available to earn interest on the investment. If v 
is low, returns will be small even if the money is invested for 
a number of years. 

Fig. 5 shows how different parameters affect the percent 
savings provided by NPC. The x-axis is the energy payment 
rate (p, which is the cost a of a Joule multiplied by the rate p 
at which the network consumes energy, and has units of $/s. 
The low and high values for node hardware and labor costs 
are discussed in Sec. VTA, while the percent savings metric 
is explained at the beginning of Sec. VI. 

NPC minimization is most effective when a large portion 
of the budget is dedicated to energy. In Fig. 5 the percent 
savings approaches 0 for values of p approaching 0 $/s. Percent 
savings are higher when /3 and (p are low. We can also see 
that the impact of changes in p, f3, and tp on the percent 
savings is significantly reduced when the MTBF w is low. NPC 
minimization works by deferring expenses to take advantage of 
interest rates, and in our scenario energy is the only deferrable 
expenditure, so the greater the portion of the budget dedicated 
to energy, the greater the percent savings. 

E. Cost Breakdown of NPC Minimization 

In this section we take a closer look as to how the cost /3 of 
node hardware, cost of labor, MTBF w, operational lifetime 
L, and interest rate v affect the percent savings and the optimal 
number of visits. We show that NPC minimization works by 
balancing decreasing energy costs with increasing labor costs, 
and study how adjusting each parameter impacts the optimal 
NPC and number of visits required to achieve this NPC. 

Figs. 6-9 show the normalized NPC for various numbers of 
visits. Each bar is segmented to show the expenditures on node 
hardware, failure repair, energy, and labor for each number of 


Fig. 6. Normalized NPC for different numbers of visits when pai'ameters 
are set to their default values. 



Fig. 7. Normalized NPC for different numbers of visits when pai'ameters 
are set to their default values, except the cost ^ of a node which is $100. 

visits being performed. Note that the labor bar segments only 
include the cost of labor for restoring energy; costs related 
to labor for failure repair are included in the failure repair 
bar segments. NPC was normalized to the cost of a network 
when K = \ for the default parameter values described at the 
end of Sec. VI-A. Fig. 6 shows the normalized NPC at the 
default parameter values. Figs. 7-9 show the normalized NPC 
where all parameters are at their default values except one. 
The adjusted value is given in the figure caption. 

Fig. 6 shows that NPC minimization balances decreasing 
energy costs with increasing labor costs. As the number 
of visits increases, the total cost of energy decreases at a 
decreasing rate while the total cost of labor appears to increase 
steadily. As nodes are only purchased once during the first 
visit, the total cost of nodes does not change. The more 
a parameter affects the rate that energy costs decrease or 
labor costs increase, the greater that parameter’s effect on the 
optimal NPC and the number of visits required to achieve it. 

Fig. 7 shows that increasing the cost of node hardware 
affects the NPC for each value of K equally. Both the 
optimal NPC and the cost when performing only one visit 
increase by the same amount. The reduction in percent savings 
from increasing /3 comes from increasing the NPC when 
K = 1 without changing the difference between it and the 
optimal NPC. The money saved on energy by performing NPC 
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Fig. 8. Normalized NPC for different numbers of visits when pai'ameters 
are set to their default values, except the cost of labor which is 140$/visit. 
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Fig. 9. Normalized NPC for different numbers of visits when parameters 
are set to their default values, except the MTBF lo which is 1 month. 


minimization remains the same as in Fig. 6. 

Fig. 8 shows that decreasing labor costs ( allows a network 
operator to perform more visits and reduce energy costs 
further, improving the percent savings. Even performing a 
large number of visits barely increases the total labor cost. 
This means that a larger number of visits can be made, and 
energy costs can be reduced further, before the total cost of 
labor exceeds the savings from NPC minimization. 

Fig. 9 shows that increasing the failure rate causes the repair 
payments to increase uniformly. The number of visits to restore 
energy does not affect network reliability, explaining why the 
repair payments increase uniformly regardless of K. A smaller 
percent of the budget is spent on energy, reducing percent 
savings without changing the optimal number of visits. 

Low interest rates v and low lifetimes L rob NPC minimiza¬ 
tion of its ability to lower energy costs. Returns on invested 
money at low interest rates and lifetimes are negligible, 
reducing the savings on future energy costs. When increasing 
K, the cost of labor immediately exceeds the savings in energy 
costs, making a single visit optimal when either i or u is low. 

VII. Conclusion 

We have addressed the problem of minimizing the Net 
Present Cost (NPC) of operating Wireless Sensor Networks 
(WSNs) by providing frameworks that determine the number 


and spacing of visits, as well as the size of their visit 
expenditures. We provided a general non-linear, non-convex 
optimization framework for minimizing NPC when the re¬ 
lationship between a visit expenditure and visit lifetime is 
known. We proposed a framework for maximizing a visit 
lifetime when given a visit expenditure, and showed that the 
relationship between visit expenditure and visit lifetime is 
linear under this framework. We developed a more efficient 
framework for minimizing NPC that takes advantage of this 
relationship, and demonstrated that equally spacing visits is a 
near-optimal strategy under such conditions. 

Compared to making a single visit, NPC minimization can 
significantly reduce costs by deferring expenditures and using 
returns on investments and network revenue to reduce future 
expenditures. It finds the optimal balance between the money 
saved on energy and the cost of labor required to deliver it to 
the network. Networks that consume relatively large amounts 
of energy, networks with long operational lifetimes, and cases 
where interest rates are high tend to benefit the most from NPC 
minimization. Large labor costs reduce the number of visits 
before the total cost of labor exceeds the savings generated 
by NPC minimization, limiting the opportunities NPC mini¬ 
mization has to lower energy costs. Low lifetimes and interest 
rates reduce the returns required by NPC minimization. 

There are a number of ways in which our NPC minimization 
framework could be applied by a network operator. In addition 
to minimizing the overall cost of a network, the network 
operator can determine the cost of deploying the network, 
the optimal number of times to visit the network to perform 
maintenance, the optimal spacing between visits, and the 
optimal visit expenditures. This information is sufficient for 
the network operator to predict the CAPEX and OPEX of the 
network. The network operator can also predict how each visit 
expenditure will be divided between node hardware, energy, 
and labor costs. A maintenance schedule that not only predicts 
when to visit the network to perform maintenance, but also 
the node hardware, energy, and labor resources required, can 
therefore be created with the NPC minimization framework. 
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